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We demonstrate that topological transport phenomena, characteristic of Weyl semimetals, namely 
the semi-quantized anomalous Hall effect and the chiral magnetic effect (equilibrium magnetic- 
field-driven current), may be thought of as two distinct manifestations of the same underlying 
phenomenon, the chiral anomaly. We show that the topological response in Weyl semimetals is 
fully described by a #-term in the action for the electromagnetic field, where 9 is not a constant 
parameter, like e.g. in topological insulators, but is a field, which has a linear dependence on the 
space-time coordinates. We also show that the #-term and the corresponding topological response 
survive for sufficiently weak translational symmetry breaking perturbations, which open a gap in 
the spectrum of the Weyl semimetal, eliminating the Weyl nodes. 



I. INTRODUCTION 

Weyl semimetals have attracted attention recently as 
a new kind of topologically-nontrivial phase of matter: 
Weyl semimetal is gapless in the bulk yet possesses pro- 
tected surface states and the corresponding topological 
transport phenomenaPEH Topological protection in this 
case results from the separation of the individual Weyl 
band-touching nodes with opposite topological charges 
in momentum space, which makes it impossible to hy- 
bridize the nodes and produce a fully gapped insulating 
state without violating translational symmetryP^^Such 
separation requires breaking of cither time-reversal (TR) 
or inversion (I) symmetry, or bothj^ as in the presence 
of TR and I all bands are doubly-degenerate by Kramers 
theorem. 

As has long been known in the quantum field the- 
ory context, chiral Weyl fermions are associated with 
the phenomenon of chiral anomalypSHHI Chiral anomaly 
manifests in nonconservation of the numbers of particles 
of a specific chirality in the presence of topologically- 
nontrivial configurations of the background gauge field 
(electromagnetic field in our context), even though these 
numbers are conserved classically (for massless particles). 
This phenomenon plays an important role in the Stan- 
dard model of particle physics I2SI22I J n the condensed mat- 
ter context, the 2+1-dimensional relative of the chiral 
anomaly, the parity anomaly, has mainly been discussed, 
due to its close relation to the quantum Hall effect EfflH 
The discovery of Weyl semimetals provides a concrete 
condensed matter system, where 3+1-dime nsional chiral 
anomaly and related effects can be realized? 11 * 13 * 14 ! 

In this paper, we focus on a specific realization of a 
Weyl semimetal in a magnetically-doped multilayer het- 
erostr ucture , made of alternating layers of topological in- 
sulatoilSSESl (TI) and normal insulator (NI) materials: 4 
This system realizes the simplest possible kind of Weyl 
semimetal, with only two Weyl nodes of opposite chiral- 
ity, the smallest number allowed by the Nielsen-Ninomiya 
theorempS in its bandstructure (identical results are ob- 
tained by magnetically doping a bulk TI with a small 
bandgap). We have demonstrated before that such a sys- 
tem possesses topologically-nontrivial transport proper- 



ties, namely the semi-quantized anomalous Hall effeclP^ 
and the chiral magnetic effect (generation of equilibrium 
current by magnetic field) P The chiral magnetic effect 
has be en kn own for some time in the particle physics 
contexlP^SHU] anc | ma y have recently been observed ex- 
perimentally in relativistic heavy ion collisions.^ Obser- 
vation of this effect in Weyl semimetals would be of sig- 
nificant interest. 

In this work we demonstrate that both the quantum 
anomalous Hall effect and the chiral magnetic effect in 
Weyl semimetals are manifestations of the same under- 
lying phenomenon, the chiral anomaly. We show that 
opposite-chirality Weyl nodes, separated in momentum 
space and in energy, give rise to an induced 0-term in the 
action of the electromagnetic field 

Se = ^J dtdr ^M)^"^^ (!) 

where h = c = 1 units are used henceforth. 0(r,t) is an 
"axion" fieldj^ which has the following form 

0(r, t) = 2b • r - 2b t, (2) 

where 2b is the separation between the Weyl nodes in 
momentum space and 2bo is the separation between the 
nodes in energy. 

In the rest of the paper, we wi ll give a derivation of 
Eq. using Fujikawa's method}- 28 * 29 * 44 * which clearly 
demonstrates the relation of the 6*-term to the chiral 
anomaly, and show that both the anomalous Hall and 
the chiral magnetic effects follow directly from Eq. 
We will also demonstrate that, somewhat contrary to the 
commonly expressed belief that Weyl semimetal is only 
topologically stable in the presence of translational sym- 
metry, which prohibits the mixing of Weyl nodes, the 
0-term in Eq. ([lj in fact survives even when the trans- 
lational symmetry is broken and the Weyl nodes are hy- 
bridized and gapped out, provided the translational sym- 
metry breaking is sufficiently weak. The quantum anoma- 
lous Hall effect and the chiral magnetic effect are thus 
more robust than the Weyl nodes themselves. 
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II. 



-TERM IN WEYL SEMIMETALS 



We start from a specific realization of a Weyl semimetal 
in a TI-NI multilayer heterostructureP The advantage of 
this system is its simplicity (and, perhaps, simplicity of 
experimental realization as well), as the Weyl semimetal 
realized in this system contains only two Weyl nodes in 
its bandstructure, i.e. the minimal number, required by 
the fermion-doubling theorem.^ The Hamiltonian, de- 
scribing the multilayer system, is given by 

H = [v f t z (z x cr) ■ kj_<5jj + ma z Sij + AgT x 5i^ 



+ -A D T + 5 jti+ i + -A d t Sj ti -i 



(3) 



The first term in Eq.([3]) describes the two (top and bot- 
tom) surface states of an individual TI layer, vp is the 
Fermi velocity, characterizing the surface Dirac fermion, 
which we take to be the same on the top and bottom sur- 
faces of each layer. kj_ is the momentum in the 2D surface 
Brillouin zone (BZ). cr is the triplet of Pauli matrices, act- 
ing on the real spin degree of freedom, and r are Pauli 
matrices, acting on the which surface pseudospin degree 
of freedom. The indices i,j label distinct TI layers. The 
second term describes exchange spin splitting of the sur- 
face states, which is induced by doping the sample with 
magnetic impurities. The remaining terms in Eq.(|3| de- 
scribe tunneling between top and bottom surfaces within 
the same TI layer (the term, proportional to As), and 
between top and bottom surfaces of neighboring TI layers 
(terms, proportional to A^). Diagonalizing Eq. ^ one 
finds, when (A s - A D ) 2 < m 2 < (A 5 + A D ) 2 , two Weyl 
nodes, separated along the z-axis in momentum space by 
a wavevector of magnitude 



2b : 



•A3, 



2A S A D 



(4) 



where d is the multilayer period. 

The multilayer model Eq. ^ possesses inversion sym- 
metry with respect to an inversion center, placed midway 
between the surfaces of any TI or NI layer. This sym- 
metry guarantees that the two Weyl nodes occur at the 
same energy as the corresponding symmetry operation 
interchanges the nodes with opposite chirality. In a real 
multilayer, this symmetry will likely not be present, and 
can in fact also be removed on purpose by making the 
top and bottom surfaces in each TI layer distinct by, e.g. 
creating a potential drop between them. In the absence 
of the inversion symmetry, the Weyl nodes of opposite 
chirality will then be shifted with respect to each other 
in energy, as well as in momentum. 6 This is achieved by 
a spin-orbit interaction term Xr y a z , which is allowed by 
the broken inversion symmetry. Adding this term to the 
multilayer Hamiltonian Eq. ([3]), leads to the energy sep- 
aration between the Weyl nodes of magnitude^ 



2&n 



A 

Apm 



(m 2 



c2 



m 2 )(m 2 



mz 



(5) 



where m 2 cl = (A s - A D ) 2 and m 2 2 = (A 5 + A^,) 2 . 

To generalize and simplify the subsequent considera- 
tions, we will move away from the microscopic model 
of TI-NI multilayer Eq. |3]) and introduce a correspond- 
ing low-energy model, obtained by expanding the micro- 
scopic multilayer Hamiltonian around the location of the 
Weyl nodes. The justification for this is that we want to 
describe universal phenomena, which depend only on the 
number and energy-momentum separation between the 
Weyl nodes of opposite chirality, but not on any other de- 
tails of the energy spectrum away from the Wey l no des. 
It can be demonstrated by an explicit calculation)^ that 
the properties of interest to us do not change if a fully 
microscopic model is used. 

We then obtain a generic low-energy model of two Weyl 
nodes of opposite chirality, separated in momentum space 
and in energy, described by the momentum-space Hamil- 
tonian 



H = T z cr ■ k + T z b + cr ■ b, 



(6) 



where we have absorbed the Fermi velocity (in general 
different in different directions) in the definition of mo- 
mentum. The operators r and cr now have a meaning, 
different from Eq. m. t now describes the node degree of 
freedom, while cr describes the conduction- valence band 
degree of freedom (nondegenerate conduction and valence 
bands touch at the Weyl nodes). Finally, we couple the 
electrons to an external electromagnetic field, and repre- 
sent the system in terms of the imaginary time action 



S = J drdr ^ [d T + ieA + b r z 
+ t z <j ■ (-iV + eA + br 2 )] ip, 



(7) 



where = (A , A) is the electromagnetic gauge poten- 
tial and tjp , ip are the 4-component spinor Grassman field 
variables. We have suppressed all explicit spinor indices 
in the Grassmann variables for brevity. We now make the 
following observation, that will play a crucial role in our 
analysis. The imaginary time action Eq. ([7]) possesses a 
chiral symmetry 



-ir'9/2 



(8) 



which expresses an apparent separate conservation of the 
number of fermions of left and right chirality. This sug- 
gests that the terms r z bo and cr ■ b in Eq. (f7]| can be 
eliminated by a gauge transformation: 



(9) 



where 8(r, r) = 2b ■ r — 2ib r and one should keep in 
mind that ip and tp' are not complex conjugates of each 
other, but are independent variables in the fermion path 
integral. The imaginary time action then becomes 



S = J drdr V> f [d T + ieA + t z ct ■ (-iV + eA)] tp, 



(10) 
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which describes two Weyl nodes of opposite chirality, ex- 
isting at the same point in momentum space and in en- 
ergy. This argument then leads one to the conclusion 
that the system of Weyl nodes, separated in energy and 
momentum, is equivalent to the system of two degenerate 
Weyl nodes and thus does not possess any special trans- 
port properties, which we know is incorrect. The missing 
link in the above naive argument is precisely the chiral 
anomaly: while the imaginary time action Eq. ^ does 
indeed possess the chiral symmetry, the gauge transfor- 
mation of Eq. |9]) changes not only the action itself, but 
also the measure of the path integral, representing the 
partition function of the system 



Z 



(11) 



where we assume, for simplicity, that the electromagnetic 
field does not fluctuate (our results do not depend on this 
assumption, as will be clear from the derivation below). 
As we will demonstrate, the change in the path inte- 
gral measure, induced by the chiral gauge transforma- 
tion Eq. ([9]), gives rise precisely to the additional 0-term 
in the action, given, after Wick rotation to real time, by 
Eq. 0. 

To derive the #-term we w ill use a simple modifica- 
tion of the Fujikawa's method]- 28 * 29 * 44 ! used originally to 
derive the chiral anomaly in the path integral language. 
We will provide all details of the derivation for readers 
which may not be familiar with the method. For a related 
application of this method to 3D TI see Ref. EH 

To begin, it is convenient to rewrite the imaginary time 
action Eq. Q in the standard "relativistic" notation. We 
introduce Dirac 7-matrices as 

7 = t x , Y l = ir y a^, 7 5 = -i 7 7 1 7 2 7 3 = r z . (12) 

Further defining 7 4 = —17°, we can rewrite Eq. Q as 



S 



d A x ip (<9 M + ieA^ + ib^ 5 ) ip, 



(13) 



where /1 = 1,...,4, tp = V^7°j ^4 = — ibo, and f d 4 x 
denotes integral over the 4-dimensional Euclidean space- 
time. We note that all matrices 7^ are antihermitian, 
while 7 5 is hermitian, and 7 s anticommutes with 7+ 



To eliminate the bu term from Eq. ( 13 ) we will apply 



an infinite sequence of infinitesimal chiral gauge transfor- 
mations 



(14) 



where 6{x) = 2b^x^ and the sign of the exponential in the 
second line above follows from the fact that 7 5 anticom- 
mutes with 7 . The variable s £ [0,1], whose differen- 
tial ds appears in the infinitesimal gauge transformation 
above, parametrizes the infinite sequence of the chiral 



gauge transformations of Eq. ( 14 1 



We need to find how the infinitesimal transforma- 



tion Eq. (14) changes the path integral measure DtpDip. 
Following Fujikawa W e consider the 3+1 dimensional 
Dirac operator, taken at stage s of the sequence of in- 
finitesimal chiral transformations 



= 7 Al [3 Al +^ + zMl- S ) 7 5 



(15) 



Since 7^ are antihermitian, is a hermitian operator. 
Suppose we can solve the eigenvalue problem for the her- 
mitian operator ]fi 



$<p„(x) = e n (j) n (x), 



(16) 



where e n are real eigenvalues and <p n (x) are 4-component 
spinor eigenfunctions (we suppress the spinor indices for 
brevity). We assume that <p n (x) wavefunctions can be 
normalized to unity (this requires assuming a finite space- 
time volume and taking the volume to infinity at the end) 



d 4 x^>l(x)(j) m (x) = S n 



(17) 



Given a complete set of eigenfunctions (p n (x), we can 
expand the Grassmann variables tp and ?p in the path 
integral as 



(18) 



where c„ and c„ are the new Grassmann variables. Anal- 
ogously, we can expand the transformed Grassmann fields 

= [1 - ids 9(x)j 5 /2] <t>n{x)c n = K{x)d n , 

n n 

i,'{x)=^{x)e- ldse{x) ^ /2 
= J2 <t>ni*W - ids6{x)f/2] = £ K(*K- 

n n 

(19) 

Defining the infinitesimal chiral transformation operator 
as 



we obtain 



ds- I d*x<f>* n {x)0(x)-f4> m {x), 



5>, 



mn^m • 



(20) 



(21) 



This immediately gives us the path integral Jacobian, 
corresponding to the chiral gauge transformation 



2\ _ lndct((7" 2 ) _ -2Trln([7) 



J = det(U~ z ) = e 

_ ids J d 4, xJ2„<t>n( x ) ( x h 5 <l>-n( x ) 



(22) 



Consider the quantity, appearing in the exponential in 
Eq. E2l 



I(x) = ^C(ar)7 5 0n(x). 



(23) 
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To understand the meaning of this quantity, we note that 
7 5 anticommutes with the Dirac operator This means 
that if 4>n ( x ) is an eigenfunction of Jj) with an eigenvalue 
e n , then j 5 (j) n (x) is an eigenfunction of lj) with eigen- 
value — e n . It follows that since the eigenvectors of a 
hermitian operator 0, corresponding to nondegenerate 
eigenvalues, are orthogonal, only zero eigenmodes con- 
tribute to J d xl(x). Then we obtain 



(24) 



where n± is the number of zero-mode eigenstates with 
positive (negative) eigenvalue of 7 5 (i.e. chirality). Thus, 
with a slight abuse of terminology, we can call I{x) a 
"local index" or "index density" of the Dirac operator 
0, in the sense that the integral of I(x) over the 3+1- 
dimensional space-time gives the analytical index of lj). 

To evaluate I(x) explicitly we use the standard method 
of heat kernel regularizationPSl The regularization is nec- 
essary, because, as written in Eq. (231, I(x) is poorly 



defined, since a finite result is obtained due to mutual 
cancellation of divergent contributions. We have 



I(x) = lim > 



= Jim ^C(x) 7 5 e-^/ M2 ^(x). (25) 

M— loo * — * 



The square of the Dirac operator in the exponential is 
given by 

2 = -d^ - (i - s )Vm + j \r,i y ]F^ 



+ z(1- S )[ 7 a*, 7 1V^7 5 



where D 



ieAp and we have used 



Dfj,, D v ] = ie{d ll A v - = ieF t 



(26) 



(27) 



Substituting this into Eq. ( 25 1 , and using the complete- 
ness relation 



E 



<t>n(. x )<l>n{y) = S(X - y), 



(28) 



we obtain 



I(x) 



lim 

M-Hx> 



lim 

M->oo 

d 4 k 



d 4 k 



tTy 5 e -ikx e -]p /M 2 el kx 



-tr/7 exp 



{ik^ + D^f {l-sfb^b^ 



M 2 



M 2 



Of, [y , Tl^ - [7*, 7"] W + D„)f 



(29) 



Mfc„ 



Rescaling momentum integration variable fc M 
and leaving only terms that survive the limit M — > oo and 



the trace operation (which must contain four 7-matrices 
and be proportional to 1/M 4 ), we finally obtain 



I(x) 



32 

e 2 
32vr ; 



, C ± U.1J -L ( 



fiu" a/3 • 



(30) 



As pointed out by Fujikawapl Eqs. (24 1 and (30) can be 



thought of as a local version of the Atiyah-Singcr index 
theoremP^l When integrated over space-time, Eq. ([30|) 



connects the analytical index of the Dirac operator IJ) 
with its topological indexP^ 

Substituting Eq. ( 30 1 back into the expression for the 



Jacobian of the infinitesimal chiral gauge transformation 
at "time" s, we obtain 



J 



-ids J d i x9(x)I(x) 



(31) 



To get the total contribution to the action from the Jaco- 
bian after the 6 A ,7 5 term has been fully eliminated from 



Eq. ( 13 ) we integrate over the variable s 



S e = i J ds I d A x9(x)I(x) 

»2 



32^ 



d 4 x9(x)e^ afs F^F a p. 



(32) 



The dependence of the imaginary time action of the sys- 
tem on has thus been fully transferred to S$, which 
describes completely the topological electromagnetic re- 
sponse of Weyl semimetal. After analytical continuation 
to real time r — > it, we obtain Eq. ([T]). Topological re- 
sponse of Weyl semimetals is thus described by an axion- 
type action, with the "axion field" 6(r, t), which depends 
linearly on the space-time coordinates. It is useful to 
compare this with the #-term in the action of the electro- 
magnetic field, characteristic of TIs. In that case 9 = n, 
which is the only nonzero value, consistent with TR sym- 
metry. The type of Weyl semimetal we are considering in 
this paper can be thought of as being obtained from a TI 
in which both TR and I symmetries have been broken.^1 
It is easy to see that broken TR allows for a nontrivial de- 
pendence of 9 on the spatial coordinates, while broken I 
allows for a nontrivial time dependence. Weyl semimetal 
can thus be thought of as being characterized by a #-term 
with the simplest nontrivial space and time dependence 
of the "axion field" 9. 

Integrating by parts and eliminating a total derivative 
term, we can rewrite Eq. ([I]) in the Chern-Simons forirP 



5» = 



8tt 2 



dtdTd IL 9e^ af> A„d a Ap. 



(33) 



Varying Eq 
obtain the fo 



( 33 1 with respect to the vector potential, we 
lowing expression for the current 



Jv 



2tt 2 



b^ vaP d a A p , n = 1,2,3, 



(34) 
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and 



3u 



2tt 2 



(35) 



It is easy to see that Eq. (34| represents the anomalous 
Hall effect Awhile Eq. ([35]) the chiral magnetic effect, 1 ^ 
i.e. generation of equilibrium current in response to an 
applied magnetic field. 

It is worth noting that the chiral magnetic effect and 
the anomalous Hall effect are closely related to the 
topological magnetoelectric effect, characterizing TR- 
invariant topological insulators.^ Indeed, we can rewrite 
Eq. (p5] as 



J = 



4tt 2 



d,0B. 



(36) 



Using j = dtP, where P is the electric polarization, i.e. 
identifying j with the polarization current (one of two 
types of currents, which may exist in the bulk of an in- 
sulator), we have 



4tt 2 



d t 9B, 



which gives 



47T 



B, 



(37) 



(38) 



in the TR-invariant case, when 9 — tt, which is precisely 
the quantized topological magnetoelectric effect In a 
Weyl semimetal sample the chiral magnetic effect can 
also be measured as charge polarization (voltage), arising 
in response to an applied external magnetic field (the 
voltage will of course not be universal in this case) . 

The anomalous Hall effect can similarly be related to 
another, equivalent, form for the topological magneto- 
electric effect. Indeed, Eq. (34 1 can be written as 



4tt 2 



V6» x E. 



(39) 



Identifying j with the magnetization current, j = VxM 
(second kind of current, possible in the bulk of an insu- 
lator), we obtain 



M 



-in 



E, 



(40) 



which is an equivalent form of the topological magneto- 
electric effect. 



between the Weyl nodes, will open up a gap and elimi- 
nate the nodes. Even more alarmingly, in the presence of 
a nonvanishing chiral chemical potential by, which shifts 
the left and right nodes in opposite directions in energy 
and thus creates perfectly nested electron and hole Fermi 
surfaces, the translational symmetry will be broken spon- 
taneously due to the formation of an excitonic condensate 
for arbitrarily weak electron-electron interactions. Thus, 
in this case, Weyl semimetal, strictly speaking, is never 
a ground state, and exists only at temperatures above 
the excitonic condensation transition temperature (which 
is most likely very low). In this section we show, that 
topological transport properties of the Weyl semimetal, 
which, as we have demonstrated in the previous section, 
are closely related to the chiral anomaly, in fact survive 
even when a spectral gap is opened due to either an ex- 
ternal potential, or spontaneously, as a result of electron- 
electron interactions, provided the gap is small enough. 

We will focus on the case of the spectral gap result- 
ing from the Coulomb interaction-driven formation of an 
excitonic condensate in the presence of nonzero chiral 
chemical potential (the final result should not depend on 
the origin of the gap). Let us first briefly demonstrate 
that nonzero chiral chemical potential bo leads to the 
spontaneous breaking of translational symmetry in the 
presence of electron-electron interactions. 

Adding electron-electron interactions to the Hamilto- 
nian Eq. ([6]), restricting ourselves to the lowest-energy de- 
grees of freedom near the nested electron and hole Fermi 
surfaces, enclosing the right (R) and left (L) Weyl nodes 
correspondingly, and eliminating the b • cr term by the 
chiral gauge transformation, discussed above, we obtain 
a BCS-like effective Hamiltonian 



H = 



J2 [( fc - 6 o)cL? c kA + (-* + 6 )4 



ki°kL 



k 

U 

V 



E 

k,k' 



„t J 



(41) 



Here U is the screened Coulomb interaction potential, 
whose approximate value can be estimated as 



U < lim 



47re 2 



1 



o q 2 + 2g(b )A7re 2 2g(b ) : 



(42) 



where g(e) = e 2 /2ir 2 is the density of states of a sin- 
gle Weyl cone and the factor of 2 in the denominator in 
Eq. (42) comes from the two Weyl nodes. Introducing 



excitonic order parameter 



V 



(43) 



III. EFFECT OF SPECTRAL GAP 

It is often stated that Weyl semimetal is a topologically 
stable phase, but only provided translational symmetry is 
preserved. Indeed, any potential (even random, but with 
a nonvanishing mean value), that can scatter electrons 



and decoupling the interaction term in Eq. (41) in the 
Hartree-Fock approximation, we obtain the standard 
BCS equation for A 



1 



Ug(b ) 



V aQh(v f T ? /2T) , (44) 



v /| 2 Ta 2 ~ 



6 



where £ c is a cutoff of the order of A$.d in Eq. (J3j) . This 
gives the critical temperature of the excitonic condensa- 
tion 



-2/Ug(b ) < 



(45) 



which can, in principle, be quite significant. 

We will now demonstrate that the spectral gap, which 
opens as a result of spontaneous, as described above, 
or due to an external potential, translational symmetry 
breaking, does not affect the induced 0-term in the ac- 
tion of the electromagnetic field and thus does not affect 
the topological response of the Weyl semimetal. A note 
of caution is in order here. The above statement is of 
course only true provided the gap is much smaller than 
the high-energy cutoff scale £ C) i.e. the energy scale at 
which the Weyl node dispersion starts significantly devi- 
ating from linearity. Once the gap becomes comparable 
to £ c , our conclusions, based on a low-energy model of 
a Weyl semimetal with linear Weyl node dispersion, can 
no longer be expected to hold. 

We introduce a fluctuating field A(r,r) in the imagi- 
nary time action of the Weyl semimetal 



drdr ip< [d T + ieA a + b r z + T z a ■ (-iV + eA 



br 2 



1 



A(r, r)r H 



;A*(r, 



1>. 



A(r, t) has the following form 

A(r,T) = A (r,T)e- 2lb " 



(46) 
(47) 



where Ag(r,r) is a slowly- varying envelope function. 
The field A(r, r) arises from Hubbard-Stratonovich (HS) 
decoupling of the electron-electron interaction term in 
Eq. (41) (we will leave out the term, quadratic in A, for 
brevity, but it is implicitly understood to be present). 
Rewriting the action in relativistic notation, we obtain 

S = i I d i x^ (lp L + iA*) ^— ^ 



(H> R + iA] 



1 + 7 5 



(48) 



where RL = 7 M [<9 M + ieA^ ± ib^(l - s)]. 

As before, we want to extract the term, proportional to 
bfj,, from the action, by a sequence of infinitesimal chiral 



gauge transformations (14 1. Following the same proce- 



dure as before, we consider the operator 



v = (U) L 



■iA* 



1 



r 



1+7 5 



(49) 



2 v 1X '2 

This operator is clearly not hermitian. We can, however, 
construct two hermitian operators 

^tjj = I" tA 2 1 a a* a*..ul /1 „\1 1 T 



r L +i ltl d^A* -A* 7 ^m(1-s) 
? R - i-fd^A - A-fb^l - s) 



1 



r 



A* A, 
(50) 



and 



r L -i ltl d^A* +A*rb^l-s) 
)\ - iy^A - A 7 ^(l - s) 



Y~ + A*A, 
(51) 

The two operators have identical nonnegative eigenval- 
ues, but different eigenfunctions 



e 2 n 4> n {x). 



(52) 



We then expand the transformed Grassmann variables 
i//(x) and 4>'(x) with respect to the complete sets <f>(x) 
and 4>{x) correspondingly and obtain the following ex- 
pression for the transformation Jacobian 



J 



_ ids J d 4 x [l(x)+I(x)]/2 



(53) 



where 



I{X) =Y,^)6{x) 1 ^ n {x) 



iim J2^n x h 



5 -V^V/M' 



<f> n {x), (54) 



and 



I(x) = J2K(^(xh 5 k(x) 



lim ^4>*(x)B(xh 



5 -VT></M' 



M—H>o 



4>„{x), (55) 



Performing exactly the same manipulations as in 
Eqs. (26 1- (30), it is then straightforward to show that 



the infinitesimal chiral transformation Jacobian has ex- 



actly the same form as in Eqs. (30), (31 ), i.e. the spectral 



gap term in Eq. (48) does not contribute to the chiral- 
transformation-induced #-term in the action. The rea- 
son for this is easy to understand. The terms that can 
contribute to I{x) and I(x) in the limit M — > 00 must 
be proportional to the product of four gamma-matrices, 
multiplied by the factor 1/M 4 , as all other terms get nul- 
lified by multiplying them with the 7 s matrix and taking 



the trace. By examining Eqs. (50) and (51) it is easy to 
convince oneself that such a term can only arise from the 
Ip R L part, as only these contain two gamma matrices. 
The HS terms in Eq. (48) explicitly break the chiral 



symmetry of the action by mixing the left and right 
fermions. This means that the HS terms themselves 
change under the chiral transformation Eq. (14). Indeed, 



the imaginary time action after the chiral transformation 
is given by 



S = Sg + i I d 4 x ip 
$ + iA* e 



-2164x4 1 7 _|_ iA a e 2ib ' lX ' 1 - — 



1>, 
(56) 



where = 7^(9^ + ieA^). Thus, in principle, while 
the dependence on b has indeed been eliminated from 
the fermionic part of the action, the dependence on 64 
(i.e. 60) remains, and could contribute to the electro- 
magnetic part of the action after fermions are integrated 
out. It is easy to see, however, that this does not hap- 
pen. Indeed, focusing on the fermionic part of the ac- 
tion in Eq. (56), and integrating out fermions, we simply 



obtain the Ginzburg-Landau action for the excitonic or- 
der parameter Ao, which must, based on symmetry and 
gauge invariance considerations, have the following gen- 
eral form 



S G l= J d 4 x [^A^Ao-HAop + ulAol 4 ] 



(57) 



where the coefficients g, r, and u depend on the gauge 
potential only through the field invariant F^F^, since 
the excitonic order parameter Ao is charge-neutral, and 
r ~ (T-T c )/T c , with T c given by Eq. ([45]). This clearly 
implies that the #-term in ( [57] ) is not renormalized by the 
excitonic part of the action. Physically this happens due 
to the fact that the gap-opening due to A is a low-energy 
phenomenon, while the 0-term contains the contribution 
of all filled states and thus can not be affected by A. 
The only way A could affect the #-term is through the 
chiral anomaly, which, as we have demonstrated above, 
also does not happen. 



IV. DISCUSSION AND CONCLUSIONS 

We will now provide a less formal explanation of the 
above results, which helps to understand the physical 
origin of the insensitivity of topological response in Weyl 
semimetals to opening up a spectral gap due to broken 
translational symmetry. Consider the low-energy Weyl 
semimetal Hamiltonian Eq. ([6]), to which we will add 
the node-mixing potential term later. We assume, for 
concreteness, that the vector b is along the i-direction, 
b = bz. We will also assume that an external orbital 
magnetic field B is applied to the system, along the i- 
direction as well, B = Bz. For clarity of the presentation 
we will separately consider two cases: b 7^ 0, 60 = and 
the general case b ^ 0, b ^ 0. In the first case we obtain 



H = T z (a x ir x + <T y TT y ) + T z a z k z + ba z 



(58) 



where 7r = — iV + eA is the kinetic momentum in mag- 
netic field. Introducing Landau level ladder operators 
a = Isi^x ~ i7Ti/)/v / 2 and a) — Isi^x + *f 9 )/v^i where 
Ib = 1/veB is the magnetic length, we obtain 

H = ~T z (a + a + a-ai) + a z (b + T z k z ), (59) 

where = 1/£b- This is easily diagonalized and we 
obtain the Landau level dispersion as 



3<\J 2to B n + (ak z + b) 



2 n > 1, 



(60) 



e(k z ) 




e(k z ) 




FIG. 1. (Color online) Landau level dispersion for b 7^ (top) 
and b = (bottom) cases. It is clear that in the b 7^ case 
there is an extra field-dependent electron density in the n = 
Landau level, given by Eq. (62 1. 



where s, a — ±, while the n — Landau level dispersions 
are given by 



— (ak z + b). 



(61) 



As obvious from Eq. ( 59 1 , the n = Landau levels are 



polarized downwards. Since the n > 1 Landau levels are 
particle-hole symmetric, they do not contribute to the 
anomalous Hall conductivity. The contribution of the 
n = levels can be deduced using Stfeda formulaPSl We 
note that the n = Landau levels give an extra field- 
dependent electron density, compared to the 6 = case 
(see Fig. [I]), given by 



Sn(B) = 



2b 1 

2n2TTe B 



2tt 2 



B. 



(62) 



The anomalous Hall conductivity can then be calculated 



as 



lim e- 



dSn(B) 



dB 



(63) 



Now suppose we turn on a weak periodic potential at 
wavevector Q z = 2b, hybridizing the Weyl nodes. In this 
case we need to fold Landau level dispersion to the re- 
duced first BZ —b<k z < b. Focusing on the n = 



8 



levels, the effect of the periodic potential will be to open 
a gap at the BZ boundary. It is clear, however, that 
the Hall conductivity remains unchanged, since Sn(B) 
remains unchanged after folding into the reduced BZ and 
gap opening. The expression for the Hall conductivity, 
Eq. ( 63 ) , can now be interpreted as conductance quan- 



tum erjh per period n/b of the Weyl node- hybridizing 
potential. 

Let us now consider the general case b 7^ 0, b ^ 0. In 
this case we can remove the term ba z from the Hamil- 



tonian by the chiral gauge transformation ( 14 ) , which 



gives rise to the corresponding 0-term in the action, and 
focus on the effect of the b^r z term. In the presence of 
the magnetic field in the z-direction and a uniform time- 
independent node- hybridizing potential A, the Hamilto- 
nian is given by 

H = ^T z (<r + a + a'a)) + T z <j z k z + r z b - At x . (64) 
v2 

Diagonalizing this we obtain the following Landau level 
dispersions 




n > 1, (65) 

while the n = Landau level dispersions are given by 

e « = aVfe-M'+A 2 , (66) 

where s,a = ±, as before. We can now calculate the 
current in response to the applied magnetic field (chiral 
magnetic effect j^S 



dk z d 



2n£ 2 B J_ A 2ir dk z 



eo- 



where de n /dk z is the z-component of the electron velocity 
in the n-th Landau level and A is a cutoff momentum, 
which we will take to infinity at the end. Since e nsa for 
n > 1 are even functions of k z , only the n = Landau 
level actually contributes to j z . Then we obtain 



2tt£| 



A &o- 
-A dk z 



c -(A)-eo_(-A)] 

(68) 



where the second line is true in the limit A/A — > 00. 
This coincides with Eq. (35). Note that, in agreement 



with our previous discussion, the large-momentum states 
are important for this effect, which makes it insensitive 
to low-energy phenomena, such as the presence of the 
gap A. 

In conclusion, in this work we have explicitly demon- 
strated that topological transport phenomena in Weyl 
semimetals are distinct manifestations of a single under- 
lying phenomenon, the chiral anomaly, and are described 
by a #-term in the action for the electromagnetic field, 
given by Eq. ([!]). We have demonstrated that the 6- 
term is insensitive to opening a gap in the spectrum of 
the Weyl semimetal due to broken (either by an exter- 
nal potential or spontaneously) translational symmetry, 
provided the gap is sufficiently small. 
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